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In response to the recent critical comment by M. Melendez and W. G. Hoover [arXiv:1206.0188v2] 
on our work [M. Campisi et al, Phys. Rev. Lett. 108, 250601 (2012)], we show that their molecular 
dynamics simulations do not disprove our theory but in fact convincingly corroborate it. 

PACS numbers: 02.70.Ns, 05.40.-a 67.85.-d 



In their comment pQ to our Letter [2], Melendez and 
Hoover claim that the Hamiltonian thermostat presented 
in [2], has "very unusual drawbacks" that "make it im- 
practical for many applications" . In support of this state- 
ment they present molecular dynamics simulations. We 
show here that, quite on the contrary, those simulations 
corroborate our theory. 

For convenience here we reproduce Figure 1 of the com- 
ment of Melendez and Hoover [T], see Fig. 1. The fig- 
ure shows the numerically computed energy probability 
distribution of one particle (blue symbols) and two par- 
ticles (black symbols) in a ID box interacting with a log- 
oscillator. In the caption the authors state: 

"The blue points correspond to a system of 
only one thermostated particle, where the re- 
sults failed to converge to the theoretical pre- 
diction during the simulation run." 

Melendez and Hoover make a mistake in calculating the 
"theoretical prediction" in the case of one particle. The 
theory predicts that the probability density function is 
the Gibbs distribution: 
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where Z(T) is the partition function, see Eq. (7) of our 
Letter. In the energy space the Gibbs distribution func- 
tion reads, as shown in statistical mechanics textbooks 
[3], and also mentioned in our Letter [2], 
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where Sls(Eg) is the density of states of the system. It is 
well known that the density of states of a system with a 
Hamiltonian consisting in the sum of n quadratic terms 
is of the form fls(Es) oc E^ 2 1 [3]. For one particle in a 
ID box, n = 1, and Qs(Es) ^ E s 1 . For two particles 
in a ID box, n = 2, and Qs{Es) is a constant, as we 
explicitly said in our Letter g]. Therefore the (blue) 
data in Fig. 1 from a simulation with one particle should 
be compared with 
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and the (black) data in Fig. 1 from a simulation with 
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FIG. 1: Figure 1 of the comment of Melendez and Hoover, 
Ref. [I] 



CO 





1 1 1 1 1 1 1 

analytic, one particle, Eq. (3) 




analytic, two particles, Eq. (4) 
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FIG. 2: Our own simulations with one and two thermostated 
particles in a ID box, compared with the correct theoretical 
predictions, Eq. (3) and (4), respectively. 
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two particles should be compared with 

e -Es/T 

P{E S ) = 735 F /T .„ (two particles) (4) 
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that is the red line in Fig. 1. 

Melendez and Hoover mistakenly compare the one- 
particle data sets with Eq. (4) , which pertains instead to 
the case of two particles only. 

Rather than evidencing a "very unusual drawback" of 
our work, the simulations of Melendez and Hoover cor- 
roborate our theory. See our own simulations in Fig. 2 
compared with the correct theoretical predictions, Eqs. 
(3) and (4) respectively. In comparing our simulations in 
Fig. 2, with Hoover and Melendez simulations in Fig. 1, 
note that the main difference is the scale of the vertical 
axis. This is because, unlike Melendez and Hoover, we 
have properly normalized the data so that the area below 
the curves is 1, as we did in our Letter [2]. 

In Fig. 2 of their comment |T|, Melendez and Hoover 
provide the results from a simulation of a ID chain of 
eighteen quartic oscillators, which are linearly coupled 
to two log-oscillators of same strength T. Actually that 
figure demonstrates a convergence of the particle tem- 
peratures toward the value given by the log-oscillators 
strength T, apparently in agreement with our theory. In 



Fig. 3 of their comment [T] , Melendez and Hoover provide 
the results from a similar simulation but for a nonequi- 
librium scenario with the two log-oscillators having dif- 
ferent strengths T± and T 2 . These simulations with linear 
chains, are neither sufficiently documented, nor conclu- 
sive. It is not possible to infer whether the simulations 
were done in a proper parameter regime and to draw any 
conclusions from them. The question whether and under 
which conditions log-oscillators may be employed to sim- 
ulate non-equilibrium situations is off-topic with respect 
to the focus of our Letter [5] , and needs further thorough 
investigations. 

To sum up, the claim of Melendez and Hoover that 
logarithmic oscillators "are not very useful in most prac- 
tical applications, whether simulations or experiments" 
has no scientific foundation. 
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